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We discuss the infrared structure of e + e~ — > 2 and 3 jets at NNLO in QCD perturbation theory and describe 
subtraction terms that render the separate parton- level contributions finite. As a first result, we find that the 
NNLO Cp contribution to the first moment of the Thrust distribution (1 - T) = -20 .4 ± 4. 



1. Introduction 

Hadronic event shapes and jet production ob- 
servables can be measured very accurately at LEP 
and future high energy electron positron collid- 
ers. By confronting these data with theoretical 
calculations, one can determine the strong cou- 
pling constant. Analyzing the different sources of 
error on these determinations, it becomes clear 
that the largest source of uncertainty is theoreti- 
cal and mainly due to the truncation of the per- 
turbation series at next-to-leading order (NLO). 

To improve this situation, the calculation 
of next-to-next-to-leading order (NNLO) correc- 
tions to jet observables is mandatory. For an n-jet 
observable, several ingredients are required; the 
two-loop n-parton matrix elements, the one-loop 
(n+l)-parton matrix elements and the tree level 
(n+2)-parton matrix elements. 

In the specific case of three-jet final states in 
e + e~ annihilation the primary process is 7* — * 
qqg, the decay of a virtual photon into a quark- 
antiquark pair accompanied by a gluon. The indi- 
vidual partonic channels that contribute through 
to NNLO are shown in Table 1. In the recent 
past, enormous progress has been made in the cal- 
culation of two-loop QCD matrix elements, which 
are now known for 7* — ► qqg |1I2| . The one- loop 
matrix elements with one additional parton [3] 
and the tree-level matrix elements with two more 
partons are also known and form part of NLO 
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programs for four-jet production j 115 j . 

2. The infrared problem 

Now that all the pieces are available, all that 
remains is to combine them in a way that pro- 
duces numerically stable results for physical ob- 
servables. To achieve this, the contributions of 
the processes shown in Table 1 are weighted by jet 
functions which select three-jet final states from 
the partonic final state momenta. At a given or- 
der, all partonic multiplicity channels contribut- 
ing to this order have to be summed. However, 
each partonic channel contains infrared singular- 
ities which, after summation, cancel among each 
other. While infrared singularities from purely 
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virtual corrections are obtained immediately after 
integration over the loop momenta, their extrac- 
tion is more involved for real emission (or mixed 
real-virtual) contributions. Here, the infrared sin- 
gularities become only explicit after integrating 
the real radiation matrix elements over the phase 
space appropriate to the jet observable under con- 
sideration. 

Exactly how to accomplish this task is 
presently an open question - see Refs. |6I7I8| and 
the problem of integrating out double real emis- 
sion contributions has so far only been addressed 
in specific calculations |9I1()I11I12] , each of which 
requires a subset of the ingredients needed for 
generic jet observables at NNLO. 

The infrared singularities of the real radiation 
contributions can be extracted using infrared sub- 
traction terms. These subtraction terms are con- 
structed such that they approximate the full real 
radiation matrix elements in all singular limits 
while still being sufficiently simple to be inte- 
grated analytically over a section of phase space 
that encompasses all regions corresponding to sin- 
gular configurations. In all cases, the subtrac- 
tion terms must be local in phase space. How- 
ever, there are two distinct approaches to de- 
rive the subtraction terms. In the first |6I7| . 
one uses phase space remappings together with 
the iterated sector decomposition |13I14| to ex- 



tract the singularities from individual terms in 
the matrix elements in terms of plus prescrip- 
tions |15ll6ll7tTH| . In the second [S], one iden- 
tifies one- and two-particle subtraction functions 
that approximate the full matrix elements in all 
of the singular limits and which are sufficiently 
simple to be integrated analytically over the un- 
resolved phase space. It is this latter approach 
that we follow here. 

One-particle subtraction at tree level is well un- 
derstood from NLO calculations and general al- 
gorithms are available for one-particle subtraction 
at one-loop |19I2UI21I22| . in a form that has re- 
cently been integrated analytically |21I22| . 

Similarly, tree-level two-particle subtraction 
terms have been extensively studied in the liter- 
ature 23 24 25 . However their integration over 
the unresolved phase space remains an outstand- 
ing issue. 



To specify the notation, we define the tree level 
m-parton contribution to the J-jet cross section 
in d-dimensions by, 



da 



d<t> m \M m \ 2 f 



(m) 



(1) 

where ~ hides all QCD-independent factors, the 
sum over all configurations with to partons and 
symmetry factors for identical partons in the fi- 
nal state. d<J> m is the phase space for m partons 
and | Al m | is the tree level m-parton matrix ele- 
ment. The jet function !Fy' defines the proce- 
dure for building J-jets out of to partons. The 
main property of J-j is that the jet observable 
defined above is collinear and infrared safe. 

3. NLO infrared subtraction terms 

At NLO, we consider the following m-jet cross 
section, 
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The cross section da^ LO has the same expression 
as the Born cross section da^ LO above except 
that to — > to+ 1, while da^ LO is the one-loop vir- 
tual correction to the m-parton Born cross section 
da B . The cross section da% LO is a local counter- 
term for da§ L0 . It has the same unintegrated 
singular behavior as da^ LO in all appropriate 
limits. Their difference is free of divergences and 
can be integrated over the (m + I )-parton phase 
space numerically. The subtraction term dcr^ LO 
has to be integrated analytically over all singu- 
lar regions of the to + 1-parton phase space. The 
resulting cross section added to the virtual con- 
tribution yields an infrared finite result. 

A systematic procedure for finding NLO in- 
frared subtraction terms in the second method is 
the dipole formalism derived by Catani and Sey- 
mour |26j . Their subtraction terms are obtained 
as sum of dipoles ^2 T>ijk (where each dipole cor- 
responds to a single infrared singular configura- 
tion) such that, 
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The dipole contribution T>ijk involves the m- 
parton amplitude depending on the redefined on- 
shell momenta pij,pk and the dipole T>ijk which 
depends only on Pi,Pj,Pk- The momenta pi, pj 
and pk are respectively the emitter, unresolved 
parton and the spectator momenta correspond- 
ing to a single dipole term. The redefined on-shcll 
momenta pij , pk are linear combinations of them. 

4. NNLO infrared subtraction terms 

At NNLO, m-jet production is induced by final 
states containing up to (m+2) partons, including 
the one- loop virtual corrections to (m + l)-parton 
final states. As at NLO, one has to introduce 
subtraction terms for the (m +1)- and (m + 2)- 
parton contributions. Schematically the NNLO 
m-jet cross section reads, 
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where da NNLO denotes the real radiation sub- 
traction term coinciding with the (m + 2)-parton 
tree level cross section da NNLO m an singular 

VS 1 

limits. Likewise, da NI j LO is the one-loop virtual 
subtraction term coinciding with the one-loop 
(m + l)-parton cross section dcrJ^ LO in all singu- 
lar limits. Finally, the two-loop correction to the 
m-parton cross section is denoted by da N ' NL0 . 

In the simple case of two-jet production, the 
subtraction terms have been fully identified |27| . 
The four-particle contribution is 
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Motivated by the fact that for three-jet produc- 
tion, the sum over dipoles is essentially equiva- 
lent to the three-parton matrix element, the four- 



parton subtraction term is given by, 
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The first term precisely cancels the real radiation 
when the four-parton configuration is perceived 
as a two-jet event. The second removes configu- 
rations when only one parton is unresolved, but 
the jet algorithm still sees only two jets. The fi- 
nal term eliminates cases when we double count. 
Taken together, Eqs. 5 and 6 yield a finite result. 
The one-loop three-parton contribution is, 
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and the subtraction term is 
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As before, we subtract the full virtual matrix el- 
ements so that the first term precisely cancels 



d<T NNLO wnen the three-parton configuration is 
perceived as a two-jet event. The second and 
third terms are merely the reappearance of the 
four-parton subtraction terms when a single par- 
ton is unresolved. Taken together, Eqs. 7 and 8 
yield a finite result. Finally, the two-parton con- 
tribution is made up of the two-loop contribution 
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and the analytically integrated subtractions 
terms (the first term on the RHS of Eqs. 6 and 
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where, the matrix elements are normalized to the 
two-parton matrix element such that 

1 



(10) 
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and 

d$3 = d$2 d$n, d<i>4 = d<I>2 d$T , (11) 

defines the dipole and tripole phase space, d<f>£) 
and d<I>T that the subtraction terms must be in- 
tegrated over. All of the master integrals neces- 
sary to perform the integrations in Eq. 10 have 
been evaluated in Ref. [TJ]. The result is an an- 
alytic expression in <i-dimensions such that when 
taken together, Eqs. 9 and 10 yield a finite re- 
sult. For the inclusive hadronic i?-ratio, = 1 
and the only remaining contribution is the two- 
parton piece which agrees with that found in the 
literature |28| . Similar results for distributions 
have been found using the sector decomposition 
method fig) . 

5. Three-jet event shapes at NNLO 

The subtraction terms based on subtracting the 
full matrix element described above are specific to 
the two-jet configurations. Subtracting the full 
tree-level five-parton matrix elements and inte- 
grating it analytically is neither sensible nor fea- 
sible. On the other hand, one might expect that 
the singular behaviour of the five-parton matrix 
elements can be represented by simpler building 
blocks that depend on only four of the five par- 
ton momenta multiplied by three-parton matrix 
elements that depend on momenta built from the 
original parton momenta. In this case, we can re- 
peat the same steps as for the two-jet subtraction 
terms and use Ref. ^B] to analytically integrate 
the subtraction term. 

At NNLO, the three-jet cross section contains 
seven colour structures, 

—&onnlo = Cf [AC\ + BCaCf + CC% (12) 
Co 

+DC A N F +EC F N F +FN F +GN Fn (-^-n\ 

Because of the QED-like behaviour of the C F 
colour factor contribution, the subtraction terms 
(for this colour factor) are very similar to the two- 
jet case. Accordingly, we have implemented these 
terms in the NLO four-jet program EERAD2 |S] 
and find that the five- and four-parton contribu- 
tions are numerically finite. At the same time, 



the analytic integration of the subtraction terms 
precisely cancels the infrared poles found in the 
two- loop matrix elements £Q. 

5.1. The average value of 1-Thrust 

One of the classic event shape variables is 
Thrust. The first moment of the Thrust distribu- 
tion (1 — T) is safe from large infrared logarithms 
and is therefore a theoretically clean and exper- 
imentally relevant observable. The perturbative 
expansion of (1 — T) is given by 

where A = 1.57, B = 32.3. The NNLO coefficient 
C receives contributions from all seven colour 
structures. We find that the Cp colour factor 
has the value, 

C*| C 2 = -20.4 ±4. 

6. Summary 

We have discussed the infrared singularity 
structure of jets in electron-positron annihilation 
at NNLO. For 2-jet events the subtraction terms 
necessary to render the individual partonic con- 
tributions finite has been worked through in de- 
tail [22]. In the case of 3-jet events, the subtrac- 
tion terms for the C F colour factor have been 
identified and analytically integrated. We have 
implemented these in a numerical program and 
produced first NNLO results for a 3-jet event 
shape. Analogous results for the remaining six 
colour structures are in progress. 
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